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ABSTRACT 

It is shown that the average magnitude - redshift diagram for a magnitude-limited 
sample of galaxies can be effectively used for mapping galaxy velocity field. This new 
approach is quantitatively introduced and applied to the CfA, ESO/LV and ZCAT 
samples of galaxies. Strong arguments are presented for infall velocities in the Virgo 
and Fornax superclusters and in the Great Wall structure. 



Key words: galaxies: dist 
scale structure of Universe 



ances and redshifts - galaxies: luminosity function - large 
galaxies: clusters: individual: Virgo, Fornax, Great Wall 



INTRODUCTION 



The most popular tool for mapping a galaxy velocity field, 
^— ^ and especially its deviation from the universal expansion, is 

to find distances using the Tully-Fisher relation (for spiral 
^— ^ galaxies) and the D n — a relation (for elliptical galaxies) - 

see Burstein (1990) and Strauss & Willick (1995) for review. 
Q\ We propose in this paper a new method of studying the 

velocity field of galaxies, which takes the observed galac- 
+T£ tic fluxes (magnitude) and redshifts as an observational 
ijclatabase. 

Q This paper presents only the basic idea of the new 
method. More details about the method and more appli- 
ed cations to the observational data will be published in a sub- 
^ sequent paper. 

The mathematical formalism used in this paper is 
closely related to that used by Hendry & Simmons (1990). 



2 THE BASIC IDEA 

If galaxies had the same luminosities (L = L ), their ob- 
served fluxes (/) would be connected with their distances 
(r) with an unique relation: 



f = L x 



(1) 



In such a case the flux (/) would be an ideal distance indi- 
cator. 

Unfortunately, this is not true due to two, widely 
known, reasons: 
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(i) the luminosities of galaxies have a very wide distribu- 
tion function, called the luminosity function and denoted by 

(ii) every galaxy sample suffers from observational selec- 
tion effects (we will consider here the most popular flux- 
limited case: / > fum). 

However, if we take into account both of these facts, 
we can introduce a very similar relation between flux and 
distance to that which holds for an ideal case (equation 1). 

Let us now consider, the expectation value of the flux for 
a given distance E( f \ r) instead of the flux (as in equation 
1). This is the so-called conditional expectation value: see 
eg. Fisz (1963). For a flux-limited sample and under the 
assumption that the luminosity function is independent of 
the space coordinates, it is expressed by the formula: 



E(f I r) 



yhere: 



L $(£) dL 

s'L ^ dL 



1 



(2) 



Llim — flim T • 

Equation (2) has a simple intuitive justification: the first fac- 
tor on the right-hand side is equal to the expectation ("av- 
erage") value of luminosity at distance r; the lower limit of 
integration in both integrals Lu m corresponds to the lumi- 
nosity cut-off induced by the flux cut-off fi tm - The integral 
in the denominator is introduced for normalization purposes. 

Equation (2) is valid for an inhomogeneous space distri- 
bution of galaxies - see APPENDIX A for the mathematical 
derivation. 

Equation (2) is equivalent to the formulae used by Segal 
(1989) and Hendry & Simmons (1990) in different contexts. 
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It is worth noting that equation (2) may by regarded 
(and used) as a new, inhomegeneity-free tool for the deter- 
mining of the luminosity function of galaxies. 

The magnitude version of the equation (2) has the form: 

f M|,m M <Z>(M)dM 

£|m|r| - c~ + 51 °" r '+ 25 131 

where: 

Mum = mum - 51og(r) - 25 

and where m and M denote the observed and absolute mag- 
nitude respectively. 



2.1 Special case: Hubble law 

In order to confront equations (2) or (3) with observations 
(namely: magnitude and redshifts for a magnitude-limited 
sample) we have to derive analogous relations as in equations 
(2) and (3) but for redshift (cz). 

In order to do this we assume that Hubble Law holds: 



cz = H r . 

We have, from equation (2): 



Eh (/ | cz) 



(4) 



(5) 



/~. $(£) dL (cz f 
where: 

Llim — flim (^^) 

and from equation (3): 

f M " m M<S>(M)dM cz 
E H (m | cz) = °! M +51og( — ) + 25 (6) 

where: 

Mum = mum — 5 log(-^-) - 25. 

H 

The subscript "H" in equations (5) and (6) indicates that 
the Hubble Law has been assumed. 



the 14.5 magnitude limit. Magnitudes in CfA are mostly in 
the Zwicky system (mz) while for ESO/LV we have choosen 
the total blue magnitude (Bt)- Raw data for both samples 
are presented in Figs 1 and 2. 

We have additionally used a subsample of the ZCAT 
redshift catalogue (Huchra et al. 1993) cut to the 15.5 mag- 
nitude limit. Since this subsample is not complete to a spec- 
ified magnitude limit we will use it in this paper in a limited 
scope. 

As observational representation of the expectation value 
of the observed magnitude m for a given redshift cz - E(m \ 
cz) the averages of the observed magnitude < m | cz > have 
been computed within 200 km/s bins in cz. 

The resultant averages are presented in Figs 3 - 7 as 
one standard deviation error bars connected by broken line. 

The same diagrams, namely: the average magnitude - 
redshift diagrams have been constructed by Segal (1989) but 
in a different theoretical "environment". 



3.2 Theory: E H (m \ cz) 

The theoretical shape of the function Eu(m | cz) has been 
computed, using equation (6), assuming the luminosity func- 
tion in the Schechter (1976) form: 



$(M) = 10°' 4 ( M *- M )( 1 + a W(- 



10' 



0.4(M» -M)\ 



(7) 



with parameters taken from Choloniewski (1982): a = —1.1 
and M t = -19.2 for the CfA sample and M t = -19.6 for 
the ESO/LV sample (H = 100). The difference between 
these two M* values reflects the systematic 0.4 magnitude 
shift between mz and Bt (Huchra et al. 1993). 

The resultant function Eu(m | cz) is insensitive to the 
assumed value of the Hubble constant H a - 

The shape of the function Eu(m \ cz) is presented in 
Figs 3 - 6 as a dotted line. We can easily see that the aver- 
age magnitude and redshift (or distance) are positively cor- 
related. 

We have not computed (nor presented) the theoretical 
curve for the ZCAT sample since this catalogue is not com- 
plete to any fixed magnitude limit. 



3 CONFRONTATION WITH OBSERVATIONS 

Let us now compare theoretical predictions expressed in 
equation (6) with the observational data. In order to do this 
we will compute the average magnitude for a given redshift: 
< m | cz >. Since the average value is always an unbiased es- 
timator of the expectation value, both functions: Eu(m \ cz) 
and < m | cz > should have, within the limits of accuracy, 
the same shape. Since Eu(ni \ cz) is computed under the 
assumption that Hubble law holds, any differences between 
them may be interpreted as a result of non-Hubble motions. 

3.1 Observations: < m | cz > 

We used as an observational database two magnitude- 
limited samples of galaxies with redshifts: the CfA sam- 
ple (Davies & Huchra 1982) for the northern sky, and a 
subsample of the ESO/LV catalogue (Lauberts & Valentijn 
1989) for the southern sky. Both samples are complete to 



3.3 Comparison and interpretation 

Both functions: theoretical - Eu(ni \ cz) and observational 
- < m | cz > are presented in Figs 3 (northern sky) and 4 
(southern sky). 

The approximate agreement between them can be seen. 
This means that the Hubble law (and other assumptions) are 
approximately correct. There are, however, some significant 
differences. 

In order to discover the source of these differences let us 
now study the same diagrams as above but for two pencil- 
beam subsamples centered on the two neighbouring super- 
clusters: Virgo (Fig. 5) and Fornax (Fig. 6). Additionally, we 
present the Great Wall region (Geller & Huchra 1989) us- 
ing data from the ZCAT catalogue (Fig. 7). For both su- 
perclusters and for the Great Wall very large differences be- 
tween theory and observations can be seen. These differences 
are correlated with the observed number density of galax- 
ies: they are largest where the observed number density is 
maximal. 
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The average magnitude and redshift are, in the region of 
the highest density of galaxies, negatively correlated. Since, 
as have been shown in subsection 3.2, the average magni- 
tude is positively correlated with distance, it suggests that 
galaxies which have redshifts larger than the average redshift 
of the supercluster have distances smaller than the average 
distance of the supercluster. And vice-versa: galaxies which 
have redshifts smaller than the average redshift of the su- 
percluster have distances larger than the average distance of 
the supercluster. 

This can be interpreted as a result of large, non-Hubble, 
infall velocities in the Virgo, Fornax and Great Wall re- 
gions. The amplitude of these velocities should be larger 
than the supercluster diameter (expressed in velocity terms) 
wich suggests that a certain fraction of galaxies inside su- 
percluters have negative relative velocities. 



4 CONCLUDING REMARKS 

We have shown that the average magnitude-redshift dia- 
gram can provide information about peculiar motion inside 
superclusters of galaxies. The same method can be also used 
for determining distances to superclusters as a whole. 

The approach introduced explicitely controls the se- 
lection effects and is insensitive to the inhomogeneities in 
space distribution. For methods which are based on the 
Tully-Fisher relation or the D n - a relation (see Burstein 
f990, and Strauss & Willick f995) these two factors influ- 
ence the results in a way which requires, in order to make 
them unbiased, quite complicated statistical treatment (see: 
Sandage f 994, Hendry & Simmons f 994, Triay, Lachieze-Rey 
& Rauzy 1994, Hudson 1994) 

The method introduced in this paper can be applied to 
many existing and future magnitude (or diameter) -limited 
samples of galaxies with redshifts. 

We encourage authors who develop models of the galaxy 
velocity field to produce also (among various others obser- 
vational predictions) the relation E(m \ cz) versus cz; which 
can be, as we show in this paper, sensitive to non-Hubble 
motions. Due to the same reason, we also encourage authors 
who publish galaxy redshifts (for a complete magnitude- 
limited samples) to construct average magnitude-redshift 
diagrams. 

One obvious modification of the method proposed in 
this paper is to use a A'araefer-limited sample instead of 
a magnitude-limited sample, and to replace galaxy mag- 
nitudes by galaxy diameters. Such an approach would be 
useful especially in a case where one wishes to study large- 
scale peculiar bulk motions of the superclusters (as a whole) 
which have large angular separations across the sky. In such 
cases the average magnitude would be strongly influenced by 
extinction in our Galaxy, while average effective (or charac- 
teristic) diameters would be completely insensitive to the 
extinction. 

We have used in our analysis the shape of the lumi- 
nosity function taken from the literature. It should be re- 
membered, however, that this shape has been computed un- 
der the assumption that the Hubble law holds. Since, as we 
have shown in this paper, there are significant deviations 
from this law, the shape of the luminosity function should 
be recomputed without assumption about linear (Hubble) 



expansion of galaxies. We will try to solve this problem in a 
future work. 
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Figure 1. The magnitude— redshift diagram for CfA sample of galaxies. 
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Figure 2. The same as Fig. 1 but for ESO/LV sample of galaxies. 
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Figure 3. The average magnitude as a function of redshift for CfA sample of galaxies. The averages are represented by one standard 
deviation error bars connected by a broken line. The dotted line represents the expected relation for the pure Hubble law. 
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Figure 4. The same as Fig. 3 but for ESO/LV sample of galaxies. 
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Figure 5. Left panel: the same as Fig. 4 but for the Virgo region (175° < a < 200°, 0° < 8 < 15°). The right panel shows the number 
of galaxies in each, 200 km/s bin used in the left panel. Note that for the Virgo cluster center (cz 1000A;m/s) the average magnitude 
and redshift are negatively correlated. This effect can be explained by infall on to this cluster. The data taken from CfA. 
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Figure 6. The same as Fig. 5 but for the Fornax region (45° < a < 60°, —40° < 8 < —30°). Note that for the Fornax cluster center 
(cz 1300 &m/s) the average magnitude and redshift are negatively correlated. The data taken from ESO/LV. 
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Figure 7. The same as Fig. 6 but for the Great Wall region (0° < a < 50°, -90° < 8 < 90°). Note that for the Great Wall center 
(cz ~ 5300A;m/s) the average magnitude and redshift are negatively correlated. The data taken from ZCAT. 
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APPENDIX A: THE DERIVATION OF 
EQUATION (2) 

Let us consider an observed joint distribution function of lu- 
minosity (L) and space coordinates (r): N(L,r) for a flux- 
limited sample of galaxies. According to the general rule 
described by Choloniewski (1991) this formula can be ex- 
pressed as: 

N(L,r) = $(i,r) S(L,r) (Al) 

where $(£, r) describes a selection-free distribution function 
of L and r and where S(L, r) denotes a selection function - 
the function which is equal to the probability that a galaxy 
having properties L and r belongs to the sample. For flux- 
limited sample we have: 

S(L, r ) = I(f>f lim ) (A2) 
where the indicator function /(.) is defined as: 

j, , I 1 if the statement is true (A3) 

I if the statement is false 

Let us assume now that: 

(i) L and r are statistically independent: 

$(i,r) = $(£) p(r) (A4) 

where $(£) denotes the luminosity function and p(r) de- 
notes the number density of galaxies; 

(ii) the space has Euclidean properties: 

/ = ( A5 ) 
Using these two assumptions in equation (Al) we have: 

N(L, r) = $(£) p(r) I(^ > f lim ). (A6) 

The next two steps are: the transformation of the variables 
from rectangular coordinates to spherical coordinates and 
the integration over angular coordinates which gives: 

N(L, r) = $(£) p(r) r 2 > f lim ) (A7) 

where p(r) denotes the average number density of galaxies 
at distance r. 

Let us now define the conditional distribution function: 

N{L 1 T) = ^ 

which gives: 

N{L ' = F*iL)I&>f« m UL { } 

and which allows us to compute the conditional expectation 
value: 

f°° T 

E(j\r)= — N(L | r)dL (AlO) 
Jo r 

and finally: 

f °° L $(£) AL i 

E(f I r) = x 4 (All) 

W 1 ' f™ §(L)dL r 2 1 ' 

where: 

Llim = flim f • (-^-1^) 



